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the coordinates. 
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Introduction 



For some stationary and non-stationary systems the time-dependent integrals of mo- 
tion have been constructed explicitly [1-5] . On the other hand it is well known that 
the integrals of the motion of classical and quantum systems are related to the sym- 
metry of the system. This relation is expressed by Noether's theorem [6,7]. The 
generalized stationary two-dimensional oscillator has been analyzed from the point 
of view of Noether's theorem [8]. In this work on the base of Noether's theorem the 
time-independent integrals of the motion, which are polynomials in the position and 
momentum operators, were found. 

Lewis and Riesenfeld [1] constructed integrals of motion quadratic in position 
and momentum operators for a time-dependent one-dimensional oscillator and for a 
charged spinless particle moving in an homogeneous varying magnetic field. In Ref. 
[2] the time-dependent integrals of motion which are linear in position and momen- 
tum operators have been found for a charged particle moving in a time-dependent 
homogeneous magnetic field. In Ref. [3] the linear time-dependent integrals of motion 
have been obtained for the quantum forced oscillator with time-dependent frequency. 
Using these integrals of motion the dynamical symmetry SU(2) and SU(1, 1) [3] 
has been associated for the charged particle moving in a varying magnetic field and 
for the parametric quantum oscillator. In Ref. [4] the time-dependent invariants 
linear in position and momentum operators have been found for multi-dimensional 
quantum systems whose hamiltonian is a non-stationary general quadratic form in 
position and momentum operators. These integrals of motion have been used to con- 
struct the propagator, coherent states and transition amplitudes between the energy 
levels of the system in terms of symplectic transformations ISp(2N,H) [4,5]. The 
time-dependent integrals of motion give an useful method to find the propagator of 
quantum systems using the system of equations found in Ref. [9] (see also [5]). Ur- 
rutia and Hernandez [10] considered a non-stationary damped harmonic oscillator, 
and have proved that the relation of linear time-dependent integrals of motion to the 
propagator is close in spirit to the Schwinger action principle method. In Ref. [11], 
and more recently [12], have been emphasized the connection of time-dependent inte- 
grals of motion of some examples of stationary physical problems with the dynamical 
symmetry concept of quantum systems. 
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On the other hand up to now it was not clear how the time-dependent integrals of 
motion found in Ref. [1-5] can be obtained from the canonical procedure of Noether's 
approach and what is the symmetry which corresponds to linear time- dependent inte- 
grals of motion of non-stationary multi-dimensional forced harmonic oscillator. The 
aim of the present work is to show that there exists such a variation of coordinates 
for the generalized oscillator [8] as well as for the non-stationary one-dimensional 
and multi-dimensional oscillators for which lagrangian variation is reduced to a to- 
tal time-derivative term. So from normal Noether's theorem we will obtain linear 
time-dependent integrals of the motion for the generalized two-dimensional oscilla- 
tor and connect these integrals with the analogue of Bragg scattering relation giving 
the time-independent invariants found in Ref. [8]. We will consider also explicitly 
the non-stationary generalized two-dimensional oscillator with varying mass and a 
time-dependent coupling term proportional to the third component of the angular 
momentum operator. It is also known that time-dependent hamiltonians [5] can 
generate squeezed states and therefore we are going to study the two-mode squeez- 
ing [13-15] in the frame of this model. Finally we show that the analysis made 
for the generalized two-dimensional oscillator can be extended for multi-dimensional 
non-stationary quadratic systems. Thus we find its corresponding time-dependent 
invariants, linear in position and momentum, by means of Noether's theorem. 

1. Generalized two-dimensional harmonic oscillator 

Following Ref. [8] we will remind in this section the properties of the integrals of 
motion and Noether's theorem for the two-dimensional generalized oscillator. We 
will start from the time-independent system 



where we are using dimensionless units, i.e., % = mo = ujq = 1. This hamiltonian 
is constructed by a two-dimensional harmonic oscillator plus M = X1P2 — X2P1, the 
projection of the angular momentum in the O z direction, with coupling constant A. 
For A = 1 the hamiltonian (1.1) describes the Landau problem of a charged particle 




(1.1) 
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moving in a constant magnetic field. Introducing the variables 

a\ = ^=(xi - ipj , cii = -^=(xi + ipi) , i=l,2, (1.2) 



and 

rj± = -j= (a\ ± ia£) , ^± = ^(aiT ia 2 ) , (1.3) 
which obey commutation relations 

[fa, Vb] = Sab] [fa, 6] = foa, ^] = , a, b = + , - , (1.4) 

we can rewrite the hamiltonian (1.1) in terms of these operators 

H = (1 + A) AT+ + (1 - A) AL , (1.5) 

where a constant term was neglected and N a , a = ±, denotes the number of quanta 
in direction a. If |A| > 1 the hamiltonian (1.1) has an energy spectrum unbounded 
from below. This spectrum, except for constant terms, is described by the formula 

E um = (1 + A)n + + (1 - A)n_ = v + Am (1.6) 

with \m\ = v , v — 2 , ...lorO and z/ = n++n_ , m = n+—n-, and n± = 0, 1, 2, . . .. 
The accidental degeneracy of the hamiltonian (1.5) was explained in [8] by means of 
the Noether's theorem and also the existence of time-independent integrals of motion. 
These constants of the motion, depending on the strength of the parameter A, are 
given by 

(1.7a, b) 
(1.7c, d) 
(1.7e,/) 

where the relative prime integers k\ and k 2 are connected with the rational number 
A by the formula 

^ = -6^ (1.8) 
1 + A k 2 K ' 

with the number e = — 1 for |A| < 1 and e = 1 for |A| > 1. The question which we 

want to answer here is how to find the time-dependent integrals of motion, linear in 







for 


\M 


>1, 






for 




< 1 , 






for 


A = 


= Tl , 
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position and momentum, for the system and how these invariants are related to the 
integrals of motion (1.7). The system with hamiltonian (1.5) is quadratic and due 
to the results in Ref. [2-5] there is a four dimensional symplectic matrix defined by 
A, which relates the linear time-dependent integrals of motion with the position and 
momentum operators: 



V 



(1.9) 



where we associate indices 1 and 2 with the labels + and — . Besides this matrix A(t) 
satisfies the first order differential equation 



dt 



(1.10) 



£ = 



(1.11) 



with the initial condition A(0) = I4. The I4 denotes the 4x4 indentity matrix and 
the matrix £ has the form 

I 2 
-I 2 

where I2 denotes de 2 x 2 identity matrix. The B is a matrix determined by the 
hamiltonian (1.5), and, in the case considered, it is 

/1 + A \ 



B = 



(1.12) 



1-A 
1+A 
V 1-A/ 

The equation (1.10) may be easily integrated and the integrals of the motion which 

are linear forms in position and momentum may be written down in the form of 

creation and annihilation operators 



A t (t) = exp{z(l + A)t} Ai(0) , 
where we choose the initial conditions 

A 1 (0) = Z+ 



A 2 (t) =exp{z(l-A)t} A 2 (0) , (1.13) 



A 2 (0) = £-- (1.14) 

Then, the invariants (1.7) are immediately obtained from the integrals of the motion 
(1.13) by calculating the operators 

A^(t)A^(t) for |A| > 1 , (1.15a) 

A^{t)Al k2 (t) for |A| < 1 , (1.156) 

A 1 (t) = £+ forA = -l, (1.15c) 

A 2 {t)=i- forA = l. (1.15d) 
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By substituting (1.13) and (1.14) into the relations (1.15a,b) and demanding that 
these integrals of motion do not depend on time, we obtain the conditions 

fci(l + A) + fc 2 (l - A) = , for |Aj > 1 , (1.16a) 
fci(l + A) -k 2 (l - A) = , for |A| < 1 . (1.166) 

These conditions remind the Bragg relation in X-ray crystallography. Thus for some 
integers k± and k 2 the dependence on time of the integrals of motion (1.15) disappears, 
in agreement with the accidental degeneracy implied by the Eq. (1.8). 



2. Linear invariants and Noether's theorem 



We have found in the previous section the linear invariants (1.13) and (1.14) simply 
guessing their form and proving that they are integrals of motion by direct checking. 
Now we will discuss how they can follow from Noether's theorem. Before doing 
that, to illustrate the procedure we consider first how the linear invariants for a one- 
dimensional parametric oscillator found in Ref. [2,3] follow from Noether's theorem. 
The lagrangian of this oscillator has the form 

L= i --^^. (2.1) 
2 2 v ' 

Following Noether's theorem procedure used in Ref. [8] let us consider the variation 
in coordinate 

5q = h(t) . (2.2) 

We have used a specific variation depending only on an arbitrary time-dependent 
function h(t), and it is straightforward to calculate the induced variation in the 
lagrangian 

SL = q (Sq)- - (u 2 q)5q = q h(t) - u 2 qh(t) . (2.3) 

This variation can be written as a total time derivative of a function, O = hq, with 
respect to time if the function h(t) satisfies the equation 



h + uj 2 (t)h = , 
6 



(2.4) 



which is identical to the equation of motion of the physical system. The integrals of 
motion, due to Noether's theorem procedure, are determined by the function O in 
the form 

FIT, 

(2.5) 



Kit) = — 5q-fl = hq-hq . 
oq 

Considering trapped particle problems, Glauber [16] pointed out that the invariant 
of Ref. [2,3] may be obtained by comparing it with the structure of a wronskian. 
Since q = p, the expression (2.5) gives linear time-dependent integrals of the motion. 
Therefore we have obtained this invariants from Noether's theorem. 

Let us now get back to the two-dimensional generalized oscillator (1.1) as another 
and more complicated example. The lagrangian for this system has the form 

L(xi , x 2 , xi , x 2 ) = ^m(xl + xl)-^m(uQ-\ 2 )(xl + xl) + \m(x 2 x 1 -x 1 x 2 ) (2.6) 

where for later convenience we take into account explicitly the mass and the frequency 
of the oscillator, however if we put m = uq = 1 we recover the case discussed in Sec. 
1. On the basis of experience with the one-dimensional problem let us consider a 
variation of coordinates x\ and x 2 of the form 



5xi = hi(t) , 5x 2 = h 2 (t) . 

The corresponding induced variation in the lagrangian (2.6) is expressed as 

5L = m{h\ + Xh 2 )xi + m(h 2 — \hi)x 2 

— m[Xh 2 + (uq — X 2 )hi]xi + m[\h\ — (u 2 — X 2 )h 2 ]x 2 . 

It may be rewritten in the following manner: 

5L = — {m(hi + \h 2 )x\ + m(h 2 — Xhi)x 2 } 

Ci/L 



(2.7) 



— m 



m 



ih; ^> ,• / 9 (Am)' , 
hi + —hi + 2Xh 2 + {J 2 - X )h\ + -h 2 

m m 

v TO ; / 2 >2\, (Am)' , 

h 2 + —h 2 - 2Xh x + {u>l - X )h 2 - ± '-h x 

m m 



X\ 



x 2 



(2.8) 



The equations of motion for the generalized oscillator with lagrangian (2.6) are 



m 



x\ = (A — oo n )xi — 2Xx 2 xi 

m 



(Am)- 



m 



-x 2 , 



,,o 2 \ ^> • m ■ (Am)- 

x 2 — (A - uJr ) )x 2 + 2Xxi x 2 ^ x\ 

mm 



(2.9) 



Comparing the last two terms of (2.8) with (2.9), we conclude again that if the 
variations of coordinates h\{t) and h 2 (i) satisfy the system of equations of motion 
of the physical system the variation of the lagrangian (2.6) takes the form of a full 
derivative 5L = dQ/dt, with the function O defined by 

O = m(hi + \h2)xi + m(h,2 — Xh\)x2 ■ (2.10) 

Due to Noether's theorem prescription we have the integrals of motion of the system 

i(t) = 7 —h 1 + 7 —h 2 -n. 

ox i ox 2 

These invariants are linear in positions and momenta 

I(t) = (pi + \mx2)hi — mh\X\ + (j> 2 — \mxi)h 2 — mli2X2 ■ (2-11) 

Thus we explained from Noether's theorem the existence of linear time-dependent 
integrals of the motion. There are four invariants because there exist four independent 
solutions for the system of equations (2.9). These different solutions are denoted by 

(k) (k) 

a superindex: h\ , h 2 , with 1 < k < 4. This set of constants of motion can be 
rewritten in the form (1.9), with the matrix A defined by the row vectors 

A« = (hf\ -m (A<*> + Xhi k) ) , m (\hP - h^)) . (2.12) 

The constants of the motion, which we denote by pio and Xio, satisfy the initial 
conditions pio(0) = pi, p 2 o(0) = p 2 , xio(0) = x\ and £20(0) = x 2 - These imply that 
the matrix A(0) = I4, i.e., 

h 1 (0) = (l,0,0,0) T , h 1 (0) = (0,-A(0), ^r,0) T , (2.13a) 

m(U) 

h 2 (0) = (0, 1, 0, 0) T , h 2 (0) = (A(0), 0, 0, ^-f , (2.13a) 

where h a (t) = (ha\ ha , ha , h a i ^) T , a = 1, 2. To find the explicit form of the 
variations h[ k \t) and h 2 k \t), we need to solve the classical equations of motion of 
the physical system. Let us introduce the change of variables 

z = hi + ih 2 , 

(2.14) 

z* = hi — ih 2 , 
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which allows to rewrite the system of differential equations for hi and h 2 as 



z + ( — - 2iX 1 z + 



K 2 "A 2 )--(Am)- 
m 



z = 



By means of the transformation 

z(t) =m" 1/2 exp|z^ oVA(t) jw(£) , 

the equation (2.15) can be simplified to the form 

w + tt 2 w = , 

with 

O 2 = J* + . 

\2m 2m 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



The Eq.(2.17), for appropriate choices of the time dependent function O, can be 
solved and then from expressions (2.14) and (2.16) we have the general solution for 
(2.9): 



hi = — -=(wexp\i [ <1t\(t) \ + w* expl — i [ (ir(A(r)l) , 
2^m\ {J J [Jo )) 

% [ dr{X{r) 
Jo 



w exp< i 



rfrA(r) 



w exp 



(2.19) 



2iy / m 

To find the linear time-dependent integrals of the motion of the generalized 
harmonic oscillator with constant parameters, we take m = ujq = 1 and A an arbitrary 
constant. The solution of equation (2.17) is 

w(t) = A exp (it) + B exp (—it) . 

From (2.19) we have the solutions for hi and h^- 

h x (t) = \A\ cos(t + Xt + (f) A ) + \B\ cos(t - At - (p B ) 
h 2 (t) = \A\ sin(t + At + 4> A ) + \B\ sin(t - At - (j> B ) 



(2.20a) 
(2.206) 



where we denoted the complex numbers A and B in polar form C = |C|exp(z$c')- 
Taking into account the initial conditions given in (2.13) we arrive to the four inde- 
pendent solutions 

hi(t) =(cosAt cost , — sin At cost , — cos At sint , sin At sint ) T , (2.21a) 
1*2 (*) =(sinAt cost , cos At cost , — sin At sint , — cos At sint ) T , (2.216) 
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Substituting these results in (2.12), we find that the matrix A can be written in 
blocks of 2 x 2 matrices 1 < k < 4 such that 

^ ^) , A fc = ^R, (2.22a, b) 

where H\= Ha = cost and ^2 = —^3 = sini. The R is a rotation matrix by an angle 
9 = Xt , i.e., 

R=( C ° s9 a S[n6 a ) . (2.23) 
y — sin t) cos V J 

Through these expressions it is immediate to get the constants of the motion, which 
can be written down in the form of annihilation operators a^o = -^(%ko + ipko)'- 

aio(t) = exp(it){ai cos(At) + a2 sin(At)} , (2.24a) 
a>2o(t) = exp(it){— a\ sin(At) + a2 cos(At)} , (2.246) 

and their correspondent hermitean conjugates. From them we built the annihilation 
operators (1.3) and found agreement with the results (1.13) of the first section. 



3. Examples with Time-Dependent Parameters 



We study systems described by the hamiltonian 

H = \ E (- + muJ ° x t) + X ^ - x ^ ■ c 3 - 1 ) 

2 \ m } 

Through Noether's theorem we get the linear time-dependent integrals of motion, and 
from these invariants, following Ref. [2-5], we can evaluate the evolution operator, 
the associated coherent states, and the correlation matrices in pq and a) a spaces. 
In these calculations for quadratic hamiltonians, the fundamental quantity is the 
symplectic matrix A, which relates linear time-dependent integrals of motion with the 
momentum and position operators. Next we are going to consider the hamiltonian 
(3.1) for several choices of the parameters m and A, and determine the analytic 
expressions for the 2x2 submatrices Afc of A. 

First, we consider an exponentially varying mass m = mo exp {'yt}, and A 
an arbitrary function of time. These parameters imply through (2.18) that O 2 = 
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u)q — 7 2 /4, so considering 7 2 < 4cu 2 , we solve Eq.(2.17), and using (2.19) together 
with the initial conditions (2.13), we determine that the matrix A of the system has 
the structure given in (2.22) but in this case the rotation angle is 9 = J* A(r) dr and 
the fik functions take the form 

, ( ^ 7 sin fit \ 
m = exp {- 7 t/2} ( cos fit +j — — J , (3.2a) 

, 2 sin fit 

^ 2 = exp {7t/2}m ^ - , (3.26) 

1 , , , sin fit 

to = exp {--ft/2} —— , (3.2c) 

mo " 

^ 4 = exp {7t/2} (^cosfit - | ■ ( 3 - 2c 

The case 7 2 > 4cUo is obtained from expressions (3.2) by replacing 

fi — > zfii 



with fii = a/7 2 /4 — Uq, and by means of the elementary relations cosifii = coshfii 
and sinzfii/zfii = sinhfii/fii. Finally, when 7 = ±2cuo, the corresponding A& ma- 
trices are obtained from (3.2) by taking the limit when fi — > 0. 

It is important to remark that if we take 7 = 0, the functions take the form 

H\ = cos u)$t , (3.3a) 

H2 = mocoosmuot (3.36) 
1 

Hz = sincuot , (3.3c) 

m uJo 

jit4 = cos u>ot , (3.3d) 

from which, if we take mo = uiq = 1, and A a constant, we recover the result given in 
the previous section. 

Sometimes it is convenient to write the quantum invariants in the form of creation 
and annihilation operators, because the eigenfunctions of the integrals of motion 
A(t) = (Ai(t), A 2 (t)) define solutions of the time-dependent Schroedinger equation 
or coherent-type states of the system: 



A\t)) M U 
11 



M ( ; t ) • (3-4) 



The matrix M is defined through the expression M = gAg 1 , with the g matrix 
given by 

*=\(f: ft), 



V2 \-iih \li 

and / = ^ti/rrioUJo is the oscillator length. It is important to emphasize that M is 
also a symplectic matrix. It is sometimes useful to write the matrix M in terms of 
the 2x2 matrices M/-, 1 < k < 4: 

Mi = ^ ( Ai - im uJoXs + A 4 H — A 2 ) , (3.6a) 

2 V m uj J 

M 2 = \ ( -Ai + zm w A 3 + A 4 + — ^— A 2 ] , (3.66) 

with M 3 = M2 and M 4 = M^. If we consider m = 1, i.e., 7 = 0, and mo = 1, the 
operators (3.4) can be written 

Ai(t) = exp{it} {ai cos(At) + a 2 sin(At)} , (3.7a) 
A 2 (t) = exp{it} {— ai sin(At) + a 2 cos(At)} , (3.76) 

where A = J * A(r) dr/^, is the average of A during the period of time t. Taking linear 
combinations as in (1.3) we arrive to the linear time-dependent invariants 

£+(£) = exp{i(u + A)tK+ , (3.8a) 
(t) = exp{z(u;o - A)t}£_ . (3.86) 

And proceeding as in (1.15) we obtain the Bragg-like conditions 

£4(^0 + A) + k 2 (uJo - A) = , for |A| > 1 , (3.9a) 
ki(uj + A) - k 2 (uj - A) = , for |A| < 1 . (3.96) 

Substituting (2.22b) into the expressions for the correlation matrices (A. 6), given 
in the Appendix A, and using that R is an orthogonal matrix, we have immediately 
that the dispersion matrices a 2 are independent on the parameter A, and are diagonal: 

2 _ hm uJo_( 1 2 t 2 

\(rr 

m u> ' 

U 

°\ q = 2moUJo ^ + Ho^o) 2 ^}^ • (3.10c) 



a PP = — o — \ 7 \2^2 + ^4 \ l 2 , (3.10a) 

a pq = ~\ \ — — + (m u )iJ4H3 jl 2 , (3.106) 



12 



Finally taking into account that for the considered system the ///-'s functions are 
given by (3.2) we get 



pp 



(t) 



exp(7t) 



h 



m uJo 



exp(-7*) 



1 + 



1 + 



7 

2fi 2 

7 2 



_7_ 
20 



sin 2 Qt '— sin 20t 



7 



sin 2 fit + -p-sin2fit 



sin 2 fit I 2 . 



(3.11a) 
(3.116) 
(3.11c) 



pqK ' 2fi 2 

For the values of the parameters 7 = 0.1 and o;o equal to 1, 1/20 and 1/30, the 
behavior of the dispersion matrices is illustrated in the Figs. 1, 2, and 3, respectively. 
In these figures, we observe that there is squeezing for the coordinates and stretching 
for the momenta. Also it is interesting to note that a pq is a negative function; 
because it is not identical to zero, there is correlation between the coordinates and the 
momenta. If we reverse the sign of 7, the roles between the dispersion for coordinates 
and momenta are interchanged, and a pq becomes positive. In particular, for the case 
7 2 < Aujq shown in Fig. 1, the dispersion matrix a pq is a negative oscillating function. 
Finally, the change of the mass as a function of time, in arbitrary units, is displayed 
in the Fig. 4. 

Now we consider a varying mass of the form 

mo , t < ; 



m(t) — < m cosh fi t 



< t < T 
T < t ; 



(3.12) 



m {O (t - T) sinhO T + coshfio^} 2 
which is substituted into (2.18). We solve the ordinary differential equation (2.17) 
for the indicated time ranges, and obtain 

{A exp (iivot) + B exp (— iuot) , t < ; 
C exp (iClt) + D exp (-iQt) , < t < T ; 
F exp (iuot) + G exp (— iuot) , T <t ; 

where O = a/c^q — ^o- By asking continuity conditions for the function w(t) and its 
derivative in t = and t = T, we find the following relations between the constants: 



F = exp(—iu> T) 
G = exp(iujoT) 



20 
A + B 
2 

A + B 



D = 



A+B A-B 
~ 20~ 



cos OT + i — sin QT ) + 

U 



- — ( cos OT + i-j^ sinOT 



O ~ \ A-B 

cos OT — i — sin flT 

u J 2 



cos OT — i sin OT 
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Substituting the constants C, D, F, and G, into the relation for w(t), it is straight- 
forward but lengthy to build the general solutions (2.19) for (2.9). The constants A 
and B are determined through the initial conditions (2.13). 

This new system will have a A matrix, which has three different functional forms 
according to the time intervals indicated in Eq. (3.12), however in all the cases the 
Afc matrices have the structure indicated in Eq. (2.22), with an R matrix identical 
to the one of the previous example : i) The functions fik for t < are obtained by 
taking 7 = into the expressions (3.3), and similarly the corresponding dispersion 
matrix is obtained through the equation (3.11). ii) For < t < T, we have to put 
the corresponding solutions (2.19) into the general relation (2.12). The resulting (ik 
functions are 



Hi 



cos Qt 



cosh fi t 

H2 = mo{fl cosh fi t sin fit + fi cos &t sinh O t}, 
sin fit 



Ms = 



A*4 



moO coshfiot 



fi 



cosh Qot cos fit — ~^ sin fit sinh fiot 



(3.13a) 
(3.136) 
(3.13c) 

(3.13d) 



Substituting these into the expressions (3.10) we arrive to the correlation matrices 
for < t < T associated to the system (3.1), with the mass parameter of the form 
(3.12): 



a pp 



hm iV0 fO 2 fin o.~ , fin , , ^ . , ~ . 

— ( u$ ~ fi 2 ^ ) ~ ^ sinh(2fi t) sin(2fit) 



ficd 2 



qq 



pq 



+ cosh(2fi t) 

— - — sech 2 (fi t) 
2moUo 

h( fig 

2\2w fi 



sin 2 (fit) + ^|cos 2 (fit) 



-^cos 2 (fit) 

fi 2 fi 2 y ' 



■ sin(2fit) - 



wofio 



°° cos 2 (fit) 



tanh(fi t) 



fi 2 Q 2 u 

For t > T the functions are given by the following cumbersome expressions: 



(3.14a) 
(3.146) 
(3.14c) 



V2 



1 {cosfiTcos[u;o(t-T)] - — sinfiTsin[w (t - T)] 
y/m/mo I 



(3.15a) 



fi m 



sinh(fioT) 



fi 

1 

UJQ 



cos (-u; (t-T) + fiV) 
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+ 


1 + — 


cos ^ 




UJQ 





cos I u>o(t — T) + OT 



H ^— - {cosh(sinh(fi ^) } 



1 + 



n 



1 - 



Wo 



UJ 



sin (-uj (t-T) + nT J 



^3 



// 4 



_1 fu>o 



%smh(O T)<i - I + 
20 1 



sin (a;o(t-T) +OT] 

(3.156) 

sinOTcos[w (t - T)] + cosOTsin[u;o(t - T)] } , (3.15c) 
sin (uj (t-T) +QT^ 



UJ 



+ 



1-° 

UJo 



sin (u (t-T) - SIT 



+ ^-{cosh(OoT) +Q (t- T) sinh(O T)} 



n 

i + — 

UJq 



1 



cu 



cos (u (t-T) -QT 



cos (u; (£-T) + fiT) 
(3.15d) 



Substituting these results into the equations (3.10) we arrive to the correlation ma- 
trices for t >T. 

To illustrate this case we make a specific choice of the parameters for the mass. 
We consider mo = ojq = 1, Oo = 0.15 and T = 10, which is displayed in Fig. 4. 
The behavior of the dispersion matrices for this selection of parameters are shown 
in Fig. 5. In spite of the change of the mass with respect to time is different that 
in the previous example (cf. Fig. 4), the general trends for the correlation matrices 
are similar. For examples, the a pp is an increasing function of time starting from its 
minimum value at t < 0, and there is squeezing for the a qq . The main difference 
appears in the correlation a pq : in this case around the axis a pq = it is an oscillating 
function for t large enough, while in the previous cases is negative or zero for any 
time. 



4. Coherent and Fock States 



Now we are going to build a general expression for the coherent-like states of the 
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studied examples. This is carried out by solving the differential equations 

A(t)<f> o (q,t) = , (4.1) 

where 

A(t) = X p p + X q q . (4.2a) 
The X p and X q are given in terms of the A& matrices by 

X p = m*/ (^i + w w A 3 ) , (4.26) 
y/2nmoUJo 

X q = y / m uJo/2h( — A 2 + A 4 ] . (4.2c) 

V m o^o / 

The solution of (4.1) yields the vacuum state of the physical system, which takes 
the form 

$ {q,t) = c(t) exp | - ^q A; 1 A g «?} , (4.3) 

where c(t) is the normalization constant. The phase of the wavefunction & (q,t) is 
chosen to guarantee that it will be a solution of the time-dependent Schroedinger 
equation. Afterwards some calculations, the final expression for the ground state 
wavefunction is 

*o(?,f) = * exp (-4 ^q- q) . (4.4) 

In all the studied examples the matrices Xk are of the form (2.22), which was used in 
the previous expression, together with the definition of the functions [i p and \x q by 
the relations 

A*p = , n } (*A»i + m uJoiJ3) , (4.5a) 
V InmoUJo 

fi q = \Jm Q u G /2Ti ( — [i 2 + /J,4 ) ■ (4.56) 

Next, we build the unitary operator 

D(a) = exp | a • A f - a* ■ , (4.6) 

which is a power series expansion of integrals of motion and so it is also an invariant. 
Now we apply (4.6) to the vacuum wavefunction (4.4), and after using a Baker- 
Campbell-Hausdorff formula [13], the Eq.(4.1), and the action of the position and 
momentum operators, we obtain the wavefunction 

$ Q (g,t) =exp|-^L + ^ ^ a-a+^-aRq^ (q,t) . (4.7) 
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This wavefunction is the general expression for the coherent-like states in the coordi- 
nate representation of the discussed physical systems. By making the substitutions 
of the appropriate expressions for the functions fi p , \i q and R into the Eq.(4.7), we 
obtain the corresponding solutions for each case. 

The wavefunction (4.7) can be expressed in terms of the multi-dimensional Her- 
mite polynomials [17] by making use of the relation 

' ni,n 2 =0 ^ ' 

Substituting this expression into (4.7) and using the form of the coherent-like states 
in the Fock-like representation, i.e., 

a l lQ! 2 2 / l a l| 2 l a 2| 



(ni,n 2 ,t|a,t) = -^=== exp | — j , (4.9) 

we can obtain the Fock-like states in the coordinate representation: 



Mp 2 I / I 



(gjnma) = $o(?,*)H n \,„ 2 > {- — Rq) . (4.10) 

The multi-dimensional Hermite polynomial can be rewritten as a product of two 
standard one-dimensional Hermite polynomials [17]: 



if i \ ( n* \( ni+n2 )/ 2 / i 

l-^* w J = Hn ivm^\ [cosdqi+sindq2] 

where we use the explicit expression of matrix R. These Fock (4.10) and coherent 
(4.7) -like states are associated to the integrals of the motion (4.2a) and therefore 
they represent squeezed and correlated states as it was shown in the previous section. 



5. General quadratic case 



In this section we show that the time-dependent invariants, linear in position and 
momentum, can be obtained through Noether's theorem procedure. Let us consider 
an arbitrary time-dependent multidimensional forced harmonic oscillator [5] 

H=^Q a B al3 (t)Q l3 + C a Q a , (5.1) 
17 



where we have defined the vector 



Q = 



/pi\ 

Pn 

which denotes the n position and the n momentum operators, and the matrices 

A B \ „ / t 



(5.2) 



B = 



C D 



C = 



G 



(5.3) 



which characterize the quadratic form for the Hamiltonian. In Eq. (5.3), A, £?, C, 
.D, stand for n x n matrices, and F, G, for n x 1 matrices. The hermiticity of the 
hamiltonian implies that the matrix B is symmetric, and this means the following 
symmetry conditions over the four constituents n x n matrices: 



A 1 = A, B l = C, D l = D . 



(5.4) 



Expanding the hamiltonian (5.1) we obtain 
1 



H 



(A af3 p a pp + 2B a pp a qp + D a pq a q[3) + F a p a + Gpqp , 



(5.5) 



where the symmetry condition (5.4) was used. 

Making a Legendre transformation and using the relation between velocities and 
momenta we get the lagrangian of the system. Following the procedure indicated 
in Sec. 2 to get the constants of the motion for this system, let us propose an 
infinitesimal variation of coordinates given by 



5q a = h a (t) 



(5.6) 



where h(t) is an arbitrary n dimensional vector depending on time. The correspond- 
ing variation induced in the lagrangian of the system can be rewritten like a total time 
derivative of a function O, if the variation of the coordinates satisfies the differential 
equation 

fL *Kp) ' + h a (A- 1 B) af3 

- (haiCA- 1 )^)' - h a {CA- x B - D) af 3 = , 



(5.7) 
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which represents the homogeneous classical equation of motion for the coordinates of 
the system [18]. For this symmetry transformation the associated conserved quanti- 
ties, according to Noether's theorem, are given by 

J = (KIvp ~ (^" 1 5) a/3 Q / 3 - A-^Fpj h a - A'lKqp + (CA-^aphaqp 
+ J dt (A~ph a Fp + (CA~ 1 ) af3 h a F f3 + h a G^j . 

(5.8) 

There are 2n invariants because the system of equations (5.7) has 2n independent 
solutions. These 2n integrals of the motion can be rewritten in the following matrix 
form 

(5.9) 



where A is a symplectic matrix in 2n dimensions, which is given in terms of the 
solutions (5.7) and the matrices characterizing the physical system 

h (2) (h^C-h^A- 1 



A(f) = 



(5.10) 



^ h (2n) (h^C - h< 2n >) A' 1 j 

The time-dependent column vector A is given by 

A k (t) = J\t (A-lh^F p + (CA-XphWFp + h™G c 



(5.11) 



In expression (5.10), the superscript denotes the different solutions for system (5.7), 
and these vector solutions are written horizontally. The initial conditions for these 
solutions are 

_ f Sj , l<i,j<n, 



hWO) 



, n + 1 < i, j 1 < 2n , 



and for their derivatives, 



"> s ' 1 .1„(0) , // - 1 < i..j< 2n . 



(5.12a) 



(5.126) 



We have proved that linear time-dependent invariants for multi-dimensional os- 
cillators can be obtained through Noether's theorem. This is achieved by considering 
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a special variation which represents a traslation along the classical trajectory of the 
quadratic system. To guarantee the existence of analytic solutions for the integrals 
of motion, it is only necessary to solve the classical equations of motion for this 
quadratic system. Because these invariants are linear in position and momentum its 
quantization is straightforward and then by means of the theory of time-dependent 
invariants, the corresponding evolution operator and other relevant quantities can be 
determined. 

Conclusions 

In this work we have found the time-dependent integrals of the motion, linear in 
position and momentum, for the generalized two-dimensional harmonic oscillator 
using the Noether's theorem. We consider in general varying mass and coupling 
strength in the Eq. (1.1). The generators of the symmetry group, for an arbitrary 
time-dependent coupling strength, were obtained as a Bragg-like condition on the 
linear invariants. Using the model hamiltonian (1.1) for specific choices of the time- 
dependent parameters we got that the four-dimensional symplectic matrix relating 
the invariants with the position and momentum operators is of the form indicated in 
Eq. (3.10), which implies that the correlation matrices are diagonal. The dispersion 
matrices obtained in the considered examples show squeezing and correlation. We 
also have constructed the corresponding solutions of the time-dependent Schroedinger 
equation. In particular, we wrote the coherent and Fock -like states in the coordi- 
nates representation. Finally, we extended the procedure to find the time-dependent 
integrals of the motion, linear in position and momentum, for the multi-dimensional 
quadratic hamiltonian using again the Noether's theorem. The time-dependent varia- 
tions that give rise to these invariants satisfy the corresponding classical homogeneous 
equations of motion for the coordinates. 

Appendix A. Correlation Matrices 

We are going to show that by means of the matrices A/-, it is immediate to evaluate 
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the 2x2 correlation matrices for the position and momentum operators [5]. Let us 
introduce the four-vector notation 

Q = {Pi,P2,x 1 ,x 2 ) = {Q a ) , l<a<4. (Al) 
Then the dispersion matrix for generalized coordinates Q a is 



1 



<f3= ^({Qa, Q(l})-(Qa)(Qf3) , (A2) 

where {Q a , Qp} denotes the anticommutator of Q a and Qp. 

From the expression for the linear invariants in terms of the matrix A, it is 
straightforward to arrive to the following expression for the correlation matrices: 

rt^) = Kim^KM , (A3) 

where cr 2 (0) is the dispersion matrix for the initial conditions, i.e., 

° (0) - - I Q _ JL _ I I ■ (A4) 
Because A is a symplectic matrix, its inverse is given by the expression A -1 = 



— SA*S and substituting this result into the expression for cr^git), we get 



a 2 (t) = -^EA'Ect^EAE . (A5) 

For the general quadratic case the matrix A has the form indicated in (2.22a), and 
after substituting it into the last expression we obtain the correlation matrices 

Vppit) = ^ m °^°( ( mo ^ )2 + A|A 4 ^ , (A6a) 

alait) = -li—^—^i + m u X\X 3 ) , (AM) 

a \^ 1 ) = ( A *Ai + {m u Q ) 2 X\X z ) . (A6c) 

Z ITIqUJq 

For application to quantum optics, sometimes is convenient to express the 2x2 
correlation matrices in terms of the creation and annihilation operators, that is 

2 If m Uo 2 1 2 \ . i 2 t A 1 \ 

°- = 2 {— *** ~ n^ a » ) + n a ™ ' ( A7a) 

2 _ 1 / m u 2 1 _2 



2 1,— ^^W^J' (A74) 

= 2 a - - ) + ■ ( A7c) 
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Figure Captions 



Figure 1. Dispersion matrices behavior in coordinates and momenta space given by 
Eq. (3.11), for the case 7 = 0.1, and mo = 1 = ujq. 

Figure 2. Dispersion matrices behavior in coordinates and momenta space given 
by Eq. (3.11) in the limit when 7 2 — > Auq. The parameters used are 7 = 0.1 and 
m = 1. 

Figure 3. Dispersion matrices behavior in coordinates and momenta space given 
by the analytic continuation of Eq. (3.11) when 7 2 > 4cu 2 . The parameters are 
7 = 0.1, mo = 1, and uq = 1/30- 

Figure 4. Time dependence of the mass. The dashed line corresponds to m(t) = 
m exp (7^), and the solid line to Eq. (3.12). The parameters used to display the 
plot are m = 1, 7 = 0.1, fi = 0.15, and T = 10. 

Figure 5. Dispersion matrices behavior in coordinates and momenta space for the 
case (3.12), with parameters O = 0.15, T = 10, and m = 1 = ui . 
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